HIGH-ORDER FINITE ELEMENTS ON PYRAMIDS 
I: APPROXIMATION SPACES. 



NILIMA NIGAM* AND JOEL PHILLIPSt 

Abstract. We present a family of high-order conforming finite elements on pyramids, which 
have polynomial traces onto the boundary of the pyramid. We first show that it is not possible to 
use approximation spaces consisting purely of polynomials. We then introduce an infinite reference 
pyramid and describe (via pullbacks) analogs of the usual Sobolev spaces on it. Using the hexahedral 
nature of this element, we are able to construct a family of approximation spaces of arbitrary order. 
The pullbacks of these onto the finite pyramid provide the requisite approximation spaces. In a 
companion article I19| . we show shape functions and degrees of freedom for these elements, with 
a view to easy computability. The discrete spaces satisfy an exact sequence property, and contain 
high-degree polynomials; this is also established in the companion paper. 

1. Introduction. High order conforming finite elements for TJ(curl) and H(div) 
spaces based on meshes composed of tetrahedra and hexahedra were first presented 
by Nedelec, [TB]- In some circumstances, it may be desirable, or even necessary, 
to use hybrid meshes which include both tetrahedral and hexahedral elements, see 
e.g. @]. If these meshes are to avoid hanging nodes then they will, in general, contain 
pyramids. A specific example arises from hybrid the application of hybrid Yee / FEM- 
FDTD schemes for Maxwell's equations, (e.g. [SUES]) where a FDTD method (using 
cubes) is used to discretize in sub-regions with uniform coefficients and rectilinear 
geometric features, and an edge finite element discretization using tetrahedra is used 
in sub-regions where the fields vary rapidly, or near complex (curvilinear or singular) 
geometric features. Pyramidal elements are used in a layer between these two meshes. 

Let n be a pyramid with a square base defined as: 

n = U = (£^C)e]R 3 U J?7 ,C>0, £<l-C, j? < l - C}- (l.l) 

It is our aim to construct high order finite elements on a pyramid. Concretely, in 
this paper we present finite element triples, (fl,U^' k (fl), £( s )' fe ), for positive integers 
k which are unisolvent conforming finite elements for i? 1 (i7), H(cm\, ft), H(div, Q), 
L 2 (fl) respectively for s = 0, 1, 2, 3. Here U^ ,k (fl) denotes the fcth order finite dimen- 
sional approximation space for the relevant Sobolev space and the sets T,^' k are the 
associated degrees of freedom. We seek finite elements with the following properties 
P1-P3: 

PI) Compatibility: Not only should the elements be conforming, but the restriction of 
each element to its triangular and quadrilateral face(s) should match that of the 
corresponding canonical tetrahedral and hexahedral finite element. Specifically, 
what should match are the traces and exterior degrees of freedom. In other 
words, these elements satisfy the correct patching conditions on inter-element 
boundaries, [12]. 

P2) Approximation: The discrete spaces U^ ,k (Q) should allow for high-order ap- 
proximation to the spaces i? 1 (fi), if (curl, f2), etc. In particular, given a positive 
integer p, it should be possible to choose k such that all polynomials of degree p 
are contained in U^' k (fl). 
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P3) Stability: The elements satisfy a commuting diagram property: 



H r {n) — — > ir -1 (curl,fl) > i3' r - 1 (div J fi) > iT -1 ^) 

n (o)| n (i)j n< 2 >| n< 3 >| ( L2 ) 

«C).*(n) — w«'*(n) -^-> w (2) ' fc (n) — w (3) ' fc (n) 



Here Dw, denote interpolation operators induced by the degrees of freedom, 
j](s),fc anc [ r j s cnosen so that the interpolation operators are well defined. 

Gradinaru and Hiptmair [T^] constructed 'Whitney' elements satisfying properties 
PI and P3 and our family of elements includes these as the lowest order case, see 
the Appendix. In the engineering literature, Zgainski et al 8, 26^ appear to have 
discovered the same first order H (curl)-conforming element independently and also 
demonstrated a second order element. In [3] the authors describe high-order finite 
elements for but not the other spaces. Graglia et al [T3] constructed ff(curl) 

and .ff(div) elements of arbitrarily high order. Similarly Sherwin {22] demonstrated 
iJ 1 -conforming elements also satisfying properties (1) and (2). These high order 
constructions provide an explicit scheme for determining nodal basis functions; none 
of them address property (3). 

The mimetic finite difference method, originally presented in 1997 |14j and further 
developed by several authors (e.g., [HI O [5] ) seeks to develop approximations on 
polyhedral meshes and hence includes pyramids as a special case. 

The major results of this paper are: the construction of conforming discrete spaces 
U ( - 0) ' k (Q),U (1) ' k (^)M {2 ' > ' k (^) and W( 3 )' fc (fi) for arbitrary (positive integer) order k. 
We show that these spaces admit convenient Hclmholtz-like decompositions, and that 
their traces on faces and edges are consistent with traces from neighbouring elements. 
Hence property PI is satisfied. 

In a companion article, [ref] we provide a description of the degrees of freedom, 
yjO),fc anc j demonstrate unisolvency. The exterior degrees of freedom agree precisely 
with those specified by neighbouring tetrahedral or hexahedral elements. Properties 
P2 and P3 are also established in the companion article, as are explicit shape func- 
tions for these high-order finite elements. To foreshadow that discussion, we will use 
the projection-based interpolation described in [9j[l0] to solve the difficult problem of 
defining the internal degrees of freedom on a pyramid. It is possible to use projection 
based interpolation for the external degrees too, and we believe that the hp frame- 
work of which it is a part will also accommodate our element. However, this is not 
our immediate objective and the external degrees described in [T7] allow for a more 
explicit exposition. 

The use of rational functions to construct finite element approximations is not 
new. In [33] Wachspress describes their use for regions in the plane bounded by 
line segments and arcs. Several previous previous constructions of finite elements 
on a pyramid use rational functions, [El [22l [3 [26]. All of these use non-polynomial 
shape functions directly on the pyramid, possibly inspired by the Duffy transformation 
mapping the cube to the pyramid. In contrast, we use a non-linear (projective) 
mapping of the reference pyramid onto an infinite pyramid where we will exploit 
hexahedral-type symmetries. This reference infinite pyramid is introduced in Section 
1.1, and we present the approximation spaces U^' k (Qoo) on this element in Section 
2. Our construction of the finite element spaces U^ ,k (fl) is detailed in Section 3. 
We verify that our approximation spaces include the first order elements of |12j in 
Section [4] We collect shape functions associated with the faces and edges on fl^ in 
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the Appendix. 

Our starting point is an observation: given polynomial data on the faces of a 
pyramid, it is not immediately obvious if these can be extended in a polynomial- 
preserving sense to the pyramid. Other authors (TTJ [2D1 12H US] have proposed using a 
macro-element that uses piecewise polynomial functions by dividing the pyramid into 
two or four tetrahedra. The difficulty in finding bases composed of functions that are 
polynomial across the whole pyramid is not a coincidence. We show in Theorem |1.1| 
that it is impossible to construct useful i? 1 (fl) pyramidal finite elements using solely 
polynomial basis functions. We do this by demonstrating an H 1 (0) function which has 
polynomial traces on the faces of the pyramid, but which does not admit a polynomial 
interpolant within the pyramid. 



Theorem 1.1. Let be the pyramid defined in (1.1 1. Consider the function 
u : Q, — > K defined by 

Then, 

1. ue H 1 ^), 

2. u has polynomial traces on the pyramid faces, 

3. There exists a k such that for all k > k , any finite element space 
satisfying property (PI ) does not contain a polynomial function that interpo- 
lates u. 

Proof. It is straightforward to verify (1). It is easy to sec u\ v —o = — £C(£ + C — 1) 
and u = on the other faces of the pyramid. This establishes (2). 

We prove (3) by contradiction. Since Q has Lipschitz boundary, we can extend u to 
a function U <E i? 1 (M 3 ) (see, for example, [1]). Choose ko so that the corresponding 
tetrahedral element interpolates U\ v= a = — + £ — 1) exactljQ Suppose that 
we could represent u = U\q by a polynomial function p, in a manner consistent 
with property (PI). Since an T^-conforming approximation must be continuous and 
traces of U on the faces will be interpolated exactly by the neighbouring tetrahedra 
and hexahedra, we must have p = U on each face of the pyramid. 

Since U = u = on four of the faces of the pyramid, we can factorise: 

K£, V, = m + C - 1)(t? + C - 1) [r«, C) + Vs(Z, V: C)] , (1-3) 
where r and s are polynomial. Further, U = —£,((£, + ( — 1) on the face rj = and so: 

p(z, o, o = ec(c + c - i)(c - iM& o = -m + c - 1), (i.4) 

which implies that (£ — l)r(£, C) = — 1- This contradicts the polynomial nature of r.U 
A similar result is presented in |24j . where it is claimed that, under the assumption 
that shape functions must be polynomial, there exists no continuously differentiable 
conforming shape functions for the pyramid which are linear / bilinear on the faces. 

1.1. The infinite reference element: pullbacks. To construct the finite el- 
ements, we shall make use of two reference elements: the finite pyramid, f2, already 



introduced in (1.1 1, and the infinite pyramid Qqo. We will typically use the symbols 



1 With the conventional understanding that the fcth order continuous element contains all poly- 
nomials of degree k, this means let ko = 3. Note that it is in fact also possible to show that we need 
rational functions for k = 1, 2 
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Fig. 1.1. Left: The infinite pyramid Qac. Right: The finite reference pyramid O 

(x, y, z) as coordinates on the infinite pyramid and (£, rj, on the finite pyramid. The 
infinite reference pyramid is defined as 

= {x = (x, y, z) e E. 3 U 00 I x, y, z > 0, x < 1, y < 1}. (1.5) 
To associate the finite and infinite pyramids, define the bijection <fi : — > O 

<f>(x, y> z) = (^r^YT^'iT^) ' #») = (o»M), (i-6) 

which is a diffcomorphism if we restrict the domain to f2oo\oo (and the range to the 
finite pyramid with its tip removed). 

Figure [Li] shows the two pyramids. The vertical faces of the infinite pyramid lie 
in the planes y = 0, x = 1, y = 1, x = 0. We denote them as Si.Ooo > S^n^ ; £3,000 > an d 
£4,0 00 respectively, and the corresponding faces on the finite pyramid S^q — ^(S^q^). 
Let -Bf2oo refer to the base face, z = 0, of the infinite pyramid and Bq the base face 
of the finite pyramid. The vertices of the finite pyramid are denoted Vi, i — 1..5, with 
V5 the point (0, 0, 1). Denote the base edges of the finite pyramid bi — S^n fl Bq and 
the other edges = VjVb- Define E as the set of all the edges of the finite pyramid. 

The infinite pyramid will serve as a tool for the construction of the function 
spaces for the elements. We thus need to understand how to map functions between 
spaces on the finite pyramid, U^' k {Vt) and the infinite pyramid, W^ s - ),fc (r2 00 ). We 
will construct the approximation spaces on the infinite pyramid to satisfy an exact 
sequence property. To have this exact sequence property preserved on the finite 
pyramid, it is necessary that the mappings between the spaces on the finite and 
infinite pyramids commute with the grad, curl and div operators. 

In the language of differential geometry, where the elements of each space can be 
considered to be proxies for 0, 1, 2 and 3-forms, the mappings should be pullbacks. 
We shall use the same notation for each pullback map - the context will never be 
ambiguous. We point the reader to [3] for an excellent treatment of the finite element 
exterior calculus. In this paper, we will switch between referring to objects as forms 
or functions, depending on the context. Formally (because we have not yet defined 



the appropriate Sobolev spaces on the infinite pyramid): 

Vu6ff l (fi) <j)*u^uo(t), (1.7a) 

VE e i?(curl, O) (j)*E = Dcf> T ■ [E o 0], (1.7b) 

VveH(div,n) <j)*v =\D<j)\D<j>- 1 -[vocf)], (1.7c) 

VqeL 2 {n) (p*q=\D<p\[qo ( f > ), (1.7d) 

/z + 1 -x\ 

where Dip is the Jacobian matrix, , qjjp I z + 1 — y . The pullback is a 

V 1/ 

bijection and the inverse pullback, (p*) 1 is equal to (p x )*. Since z > 0, Dp T Dp is 
positive definite. 

1.2. The infinite reference element: Sobolev spaces. The infinite reference 
pyramid has obvious symmetries, which make it easier to specify and analyze approx- 
imation spaces. However, it has semi-infinite extent along the z-direction, and we 
must therefore describe analogues of 7? 1 (f2), if (curl, f2) etc. on f^. Not surprisingly, 
these Sobolev spaces will have weighted norms. 

Definition 1.2. LetQoo be the infinite pyramid defined in (1.5), andp : 0^ — > Q 
be the pullback map. We define the inner product spaces 

• H^ u (fl 00 ) is the set of scalar-valued functions v : floo — > K with inner product 



(«, «)Hi(n„) := / Ti ~ . 4 + (Vu) 1 AVvdx. 
Jn^ K L + z ) 

Here A — \Dp\Dpr 1 Dp~ lT is positive definite. 

• H w (cwl, fioo) is the set of vector-valued functions (1-forms) F : 0,^ — ¥ (R) c 
with inner product 

(F, G) Hmicnll , noo ) := / (FfA(G) + (curlF) T £(curlG)dx. 

Here B = \Dp~ 1 \Dpi T Dp, and is positive definite. 

• H w (diy,fl 00 ) is the set of vector-valued functions (2-forms) F : floo — > (M.y 
with inner product 

(F,G) Hw{div ^ } := jf (F) T 6(G) + (divF) T (l + z) 4 (divG)rfx. 

• L^lflca) is the set of scalar-valued functions (3- forms) with inner product, 

( u > v )Ll(n x ) ■= / (1 + z) 4 (uv)dx. 



We observe that the inner products in these Sobolev spaces on the infinite pyramid 
are weighted by powers of rjfpzV ^ ne subscript w is to make clear that these are 
weighted norms. It is clear that these inner product spaces are, in fact, Hilbert 
spaces. We can now identify these Hilbert spaces with the usual Sobolev spaces on 
the finite pyramid. 

Lemma 1.3. The spaces H^ l> (n oo ),H w (c\\il,n oo ),H w (div,Q 00 ) and L^ft^,) de- 
fined above are Hilbert spaces. Morever, p* : ff 1 (f2) — > H^^oo) is an isometry. The 
analogous statements are true for H w (cur\, r2 00 ),_ff 11I (div, floo) and L^fioo). 
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Proof. The pullbacks, (f>* are formally bijections because f2 and have the same 
dimension. Suppose u is a 0-form in 7^ 1 (SI) and let u = <p*u. Then 

\Hh m = f |^|Kx)| 2 rfx= / —^—\ u (x)\ 2 dx. 

Now, the gradient and pull-back operator commute. Since Vu is a one-form, we use 
the appropriate pull-back to obtain 



/ \D4>- 17 

Jn 

[ \D<j)\\D(j)- lT \7u\ 2 dx= [ \7u T A\7udx. 



Hence ||u||^i (n) = IN|| 2(n) + ||Vu|| 2 2(0) = ||u|| H i (noo ). The proofs for H w (curl, O^), 
iJ„,(div, f^oo) and L^^oo) follow analogously. □ 

The weights are entirely specified by the projective mapping, <f>, and the asso- 
ciated pull-backs for the various forms. It is important to note, for example, that 

IMl!s,(fioo) = fnoo (T^F dx if u is a zero form ' while H^lli^no.) = Jn x " 2 ( x + z ) 4dx 
if ii is a 3-form. 

We also collect here concrete instantiations of the inverse pullback mapping. 

Vuefli(fioo), (f 1 )'"^"^. (l-8a) 

/l \ 

Vf?eiT U) curl,fi 00 ), (cj>- 1 )*E=[(l + z) 1 •£]o^ 1 , (1.8b) 




V«e J ff w (div,0 oo ), = [(! + *) l + z -y -vjo^- 1 , (1.8c) 

V o o i y 

VgeL^floo), (O*9=[(l + «) 4 9]°0 _1 - (l-8d) 

1.3. Rotations and traces. Define -Rq^ : ^oo — > ^oo to be the affine mapping 
that sends the infinite pyramid to itself and rotates it a quarter turn about the axis 
x = y = \, that is, the vertical face Si.n^ is mapped to £2,0^, the face S^o^ is 
mapped to S^n^ , etc. Explicitly, 

Rn 00 --{x,y,z)^-{l-y,x,z). (1.9) 

We can also define a mapping that sends the finite pyramid to itself, rotating the 
faces, R : O — >• Cl by 

R = <j>oR noo o<j>-\ R-.^^o^a-v-ctc)- 

It is clear that if an approximation space W( s )' fc (f2oo) is invariant under the mapping 
Rnoo, its (inverse) pullback to the finite pyramid will be invariant under R. This 
property will prove convenient when we consider exterior shape functions and exterior 
degrees of freedom. 

From differential geometry, we expect that zero trace data will be preserved by 
the pullback mapping. The following lemma makes this explicit in our concrete vector 
calculus formulation. We suppose that Sn^ is a surface of the infinite pyramid and 
let Sq be its image under <f> on the finite pyramid (although the lemma is true in 
general for any pair of surfaces, S\ and S2 and bijection ip : Si — > S2). 

Lemma 1.4. 



• A 1-form u is normal to Sq at a point £ = <j){x) if and only if the pullback 
(f>*u is normal to Sq^ at x. 

• A 2-form u is tangent to Sq at a point £ = 4>(x) if and only if the pullback 
4>*u is tangent to Sq^ at x. 

Proof. Let Sq be described (locally) by Sq — {£ : /(£) = 0}. Define g = f o (p, 
then Sq^ = {x : g(x) = 0}. To establish the first result, let u be a 1-form which is 
normal to Sq at £, then 

«(0 = A(C)v/(0 (f .10) 



for some scalar function A. By the chain rule, and substituting (1.10) 

Vg(x) = (D^) T {x) ■ (V/)(*(aO) = (Dc^f(x) uWx)) " " 



\(cj>(x)) \(4>(x)) 
A(0(x))V 5 (x) = 4>*u(x). 



Hence, <p*u is normal to S at x. Since </> is a bijection, the converse is also true. To 
establish the second result, let u be a 2-form which is tangent to Sq then u ■ V/ = 0. 
The chain rule gives us Vg — (D<fi) T (V /) o (j> and by definition of the pullback, 
<f)*u = |L>0|(_D0) _1 • (u o 0), hence: 

Pu ■ V 5 = \D<f>\(u o 0) T • (i?^ 1 ) 7, • (Dcjif ■ [(V/) o 4>] 
= |£0|( U T -V/)o^O. 

Hence is tangent to Sq^. Since <p is a bijection, the converse is also true.D 

In order to satisfy property PI we need to understand the traces of the functions 
in the approximation spaces of the tetrahedral and hexahedral elements to which the 
pyramid will be patched. We will use the following notation for trace maps: 
Definition 1.5. For s = 0, 1, 2 define: 

• P|q : the trace to each vertical face, Si n^ of the infinite pyramid, for 
i = 0,1,2,3. 

• rf q.' the trace to each face Si_n of the finite pyramid, Q for i = 0, 1, 2, 3. 

• fioo an d -Tfj q. t/ie traces to Bq^ and Bq. 

When s = 0, the trace of a suitably smooth function is its point value on the 
boundary; when s = 1, it is the tangential components and when s = 2 it is the 
normal component. It is helpful to recognise that if u is a differential form defined 
on a manifold M, its trace map, on some sub-manifold, N, is the pullback of 
the inclusion map, N M. We refer the reader to ([2], pg 41 ff.) for a more 
detailed exposition. The consequence for us is that trace maps commute with <f>* , 
(e.g. Tf n o </)* — cf>* o r| n ) so results we establish on faces and edges of fioo will 
carry over to the finite pyramid. 

Using Q~7] as our source, we catalogue here the traces of the tetrahedral and 
hexahedral approximation spaces that must be matched by our spaces, U^' k (fl). 
Since our construction will start on the infinite pyramid, we also show some of the 
pullbacks of these spaces to the infinite pyramid. 

• Compatibility of traces r° it: The approximation space for the tetrahe- 
dral if ^conforming element is P k , therefore if u € U^' k (£i), it must satisfy 
T\ n u G P k [£, (]. That is, T° in u is a sum of terms of the form £ a ( c (l-() k - a ~ c 
with a + c < k. Pulled back to the infinite pyramid, we see that these are 
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of the form T° n <p*u — m3jfc € P^x, z] 2 By symmetry we need similar 



(i+^) fc 

conditions on the other vertical faces of the pyramid. The approximation 
space for the hexahedral element is Q k ' k > k and so we need T° B n u 6 Q k ' k [x, y] 
Compatibility of traces r^it: For 1— forms, the approximation space for 
the fcth order P(curl)-conforming tetrahedral element presented in [T7] is 
(P k - 1 ) 3 ®S k ' 3 , where S k ' 3 = {u e {P k ) 3 | (£-£ ) ' u = °l for some arbitrary 
base point, £ . We can take £ = (0,0, 1), and observe that the face trace 
onto Si t n must be contained in (P k ~ 1 ) 2 © S k ' 2 where S k > 2 = {w e (P fe ) 2 |(£- 
£ ) • w = 0}. This specifies the compatibility constraint the trace T\ qU must 

obey. The pullback of this space is [P k ~l[x,z]) 2 ® P^Ix, 1 + z] ( _ 

We can use symmetry to define similar conditions on the other vertical faces 
of the pyramid. The traces of fcth order hexahedral conforming elements for 
-ff(curl) lie in the polynomial space Q k ~ l k (x,y) x Q k ' k ~ 1 (x,y). 
Compatibility of traces T 2 n u: For 2— forms, F 2 qU are the normal traces of 
u onto the face Si of the pyramid. The polynomial space used for the iJ(div)- 
conforming element on the tetrahedron is (pk -1 ) 3 ^_P fc_1 . The functions 
derived from the homogeneous polynomial space make no contribution the 
face traces. Therefore, the normal traces onto the faces of a tetrahedron will 
lie in (P fc_1 ), and this is the compatibility condition we must ensure for T 2 n . 
Consider a 2— form, u, whose normal trace is U ■ n € P k ^ 1 on face S\. The 
pull-back of the associated normal vector field to the infinite pyramid will 
have the form 



I / o 

" = , — I a°2 c (l + z ) k - 1 - a - c | . a r <■ ; /V - I . 



(1 + z 



_ So ) r^f U 6P^ 2 1 [ 3: ,4 

The discussions above suggest the face-wise constraints which must be satisfied by 
any approximation spaces bt s (Q). However, as was demonstrated by Theorem 1 1 . 1 1 the 
difficulty of pyramidal interpolation stems from the need to find intcrpolants that 
match trace data on all the faces simultaneously. This point will be discussed later. 

2. The approximation spaces W^ s - I,fc (r2 00 ) on the infinite pyramid. We are 

now ready to construct the approximation spaces on the infinite pyramid. These will 
be used, via the pullback map, to construct the approximation spaces on the finite 
pyramid. As a preliminary step, we identify families of "rational polynomials" on fioo 
which will be used extensively. We want the spaces on the finite pyramid £1 to contain 
all polynomials up to a specified degree. Consider the effect of the pullback mapping 
cj> on a polynomial of degree k, p = ^"ry^C 7 € P 1 (f2), where a + f3 + 7 = k: 



From Lemma 1.3 <f>*p € if^,(f2oo). This motivates our next definition: 

Definition 2.1. Let Q l ' m,n (x,y, z) to be the space of polynomials of maximum 
degree l,m,n in x,y,z respectively. Define the space of k-weighted tensor product 



2 the space is defined in | |2,3[ ) 
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polynomials 

Q'.™."( W ) = j^-^ : u G Q<' m '"( W )}. 
It will be helpful to remember the inclusion: 

Q l,m,n c QlM\ (2 . 2) 

Let P n (x,y, z) be polynomials of maximum total degree n in (x,y,z) and define the 
space of k-weighted polynomials of degree n 

P£{x,y,z) = { U ^ z *l : ufi,^)^), (2.3) 



2.1. i7 ) J ; (ri 00 )-conforming approximation spaces. We begin our construc- 
tion by noting that in [17] . the finite element approximation space for a hexahedral 
element consists of polynomials of the form p = ^rj^C, 1 . From (2.1|, we know that 

4>*p — x ^ z ^k G i? 1 (ri 00 ), ifa + /3 + 7 = fc. We might therefore expect to base an 
approximation space for H 1 ^^,) on the fc-weighted space, Q\' k ' k . However, there 
are some elements of Q k ' kk which, when pulled back to the finite pyramid, become 

a b k 

undefined at £ = (0,0, 1). The problem arises with elements of the form ^^Jjk on 
the infinite pyramid. The following examples are illustrative. 

Example 2.2. Consider the monomial pi(x,y, z) = x on the infinite pyra- 
mid. The inverse pull-back onto the finite pyramid is ((j)~ 1 )*p = fz^- The limit 
\im^^ o ((/)~ 1 )*p depends on the path by which we approach £ . Specifically, if we take 
the path a\(t) = (A(l — £) , 0, t) then lim t _ i .i(</>~ 1 )*p(aA(£)) = A. 

Example 2.3. Consider the function P2(x, y, z) = (\+ z )k on the infinite pyramid. 
Pulled back to the finite pyramid, (4>~ 1 )*p2 — C, k ■ We must therefore retain p2 in the 
approximation space on the infinite pyramid. 

Lemma 2.4. Let f2oo be the infinite pyramid described above, and k > 1 be a fixed 
integer. 

• Functions p(x,y,z) := € Q k k ' k ' k ^ 1 satisfy p € H^Qoo). 

• If p{x,V,z) = r -ji0jh ,r(x,y) G Q k ^ k (x,y), then lim^^ -1 )*^) is only 
well-defined ifr(x,y) = 1. 

Proof. We can verify the first statement by using Definition |1.2| The second 
statement can be proved by contradiction, as in Example |2.2| □ 

This result and the examples suggest the basis functions to include in a finite- 
dimensional approximation space for H^^oo). 

Definition 2.5. Let k be a positive integer. We define the underlying spaces: 

z k 



U^ k (^oo) = Qr ^©spanj-^j. (2.4) 

Lemma 2.6. The rational polynomials | > Q < a, b < k,0 < c < k ~ l| and 

(!^. z )fc form a basis for I4(°)> k (floe) . Moreover, W( ^ fe (f2 oo ) can be represented as 

ZTM^fi^) = { u e Q k k k ' k : Vu e q^ 1 ^- 1 x qM- 1 -*- 1 x Q^ 1 }. (2.5) 



Proof. The basis functions are determined by using the definition of U^' k (fl 



and Lemma 



2.4 



The gradients of rational functions of the form jj^^k are 1-forms in 



Q k-i,k,k-i x Qfc.fc-x.fc-i x qW-i, Moreover, = (0,0, (^r) T The re- 

verse inclusion follows readily by a similar calculation. This establishes the alternative 
characterization of W(°)' fe (r2 00 ). □ 

We must now constrain these spaces to obtain the approximation spaces which 



satisfy the compatibility constraints PI. This follows the discussion in Section 1.3 

Definition 2.7. Let k be a positive integer. We define the k-th order approxi 
mation spaces W^ - ) ' fc (r2 oo ): 



W^^n*,) = {ue WW.*^) | I\ noo G P£[x,z], similarly on S ijQoa ,i = 2,3,4}. 

(2.6) 

Since we will be working in the projection-based interpolation framework while 
specifying internal degrees of freedom, we define a subspace U^' k (f^oo), consist- 
ing of functions in I4^ ,k (iloo) with zero trace on the boundary of Sloo. Clearly, 
W^'^fioo) = {x(l - x)y{\ ~ y)zu, u G Q^-^}. 

In the Appendix, we present the shape functions in W^ - ),fe (r2 oo ) associated with 
the faces, edges and vertices of ftoo. These are linearly independent. Moreover, the 
number of these functions associated with a given face is exactly the same as the 
dimension of trace spaces from neighbouring simpliccs. 

2.2. flu, (curl, r2 00 )-conforming approximation spaces. We now present the 
construction of the approximation space U'' 1 '' k (Q 0o ) of fl„,(curl, fioo)- As before, 
this construction is motivated by the ultimate goal of constructing a finite element 
approximation space for fl(curl, f2) which satisfies property (PI). 

To satisfy the commuting diagram property we will need, at the very least, to have 
VW (0),fe (fi oo ) C L/W'^fioo). The alternate characterization ofU^(Q 00 ) in Lemma 
2.6 suggests that we might consider the space Qk- 1 ^^- 1 x Q fc ' _1 ' _1 x Q k ' + {~ 1 as a 
candidate for an approximation space for fl(curl, f2oo). However, this space includes 
functions that are undefined at the point £ = (0, 0, 1) on the finite pyramid. We 
must be careful here to identify what kind of discontinuities we wish to exclude on 
the finite pyramid. We are not interested in point values of our functions, only the 
tangential components. 

/»/(! + *)\ 

Example 2.8. Consideru = G Q^- 1 ^^- 1 x q^-i^-i x qJM,*-^ 

V o J 

/ 17/(1-0 

Its (inverse) pullback to the finite pyramid is, (<fi 1 )*u = 

W(i-C) 2 , 

Let v = (0, —1, 1) and consider the path a\(t) = (A(l — t), 1 — t, t). This path lies 
on the face S3 for A G [0, 1], and S3 is tangent to v. The limit of the component of 
(<p^ 1 )*u tangent to v at £ along the path a\ is lim t _>.i u(a\(t)) ■ v = A. 

/ r x z\ 

Example 2.9. Let u = n+zWfi \ r v z ) ' r e Q k ' k [ x >y]> be a 1-form defined on 




the infinite pyramid. Note that we can write u = V( p^fp+i ) 



( \ 



(k+l)rz k - 1 



, from 



which it is apparent that u G H w (cur\, fioo)- 

With these examples in hand, we are able to define approximation spaces for 
if „, (curl, fioo). 

Definition 2.10. Let k > 1 be an integer. We define the underlying space 



pilvl, reg^i,n. (2.7) 



v fe-l /fx^ 1 

ife.fer 



.4n equivalent characterization of the underlying space U ( - 1 ^ k (fl 00 is given as 



UW^(n oo ) = {ueQ k k - i 1 k - k xQ 



fe+i ~ ^fe+i x ^fe+i 

x u G Q fc+2 x g fe+2 ' ^ X Q fe+2 ' ^ }, (2.8) 

We add surface constraints to get the full definition of the approximation space: 
let Hi be the (outward) normal to the vertical faces S^q^ of Q^. Then Tj n u := 

u x ^k.sw 

Definition 2.11. Let k > 1 fean integer. Define 

w'^fioo) - g ^ | g (p^Mx, z]) 2 © i + *] ^ +f) 

and similarly on S^.q^ , i = 2, 3, 4, }, (2.9) 

w/iere 

i + *] = (T^ny s P an M 1 + *) fc_1 ~ a > o < a < fc - 1} . 

We can denote the subset U^' k (floo) consisting of elements whose (tangential) traces 
vanish on as (Qoo). 

As for ti^ ,k (floo)i in the Appendix we have tabulated the edge and face shape 
functions for D,^. These are linearly independent, and are consistent along shared 
edges. The same will be true of the pull-backs onto the finite pyramid. 

2.3. H w (div,fl ca ) and L 2 u (r2 00 )-conforming approximation spaces. Follow- 
ing a similar strategy to the previous sections, we wish to construct approximation 
spaces U^' k (floo) for H w (div,Q oc ), such that their pull-backs to the finite pyramid 
provide approximation spaces for i?(div, 0). Again, we want curlu G W( 2 )' fe (floo), Vu G 
W( 1 )' fe (f2 (X) ). Now, the curls of functions u G WW'^floo) satisfy 

V x u G gt fc 2 _1,fc_1 x ^;2' fe,fe " 1 x Q^i'*" 1 '*- 1 . 



Not all of these will have well-defined normal traces, and we must exclude these. 
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Definition 2.12. The underlying space for the H(div)- conforming element is 
defined as: 



\ fc-i / 2* 
2s © 



(2.10) 



Here s(x,y) £ Q k 1,k [x,y], t(x,y) € Q k,k l [x,y]. An alternate characterization of 
this space is 

u^^oo) = {u e Ql , +2 1 ' k ~ 1 x Qi+l' k ' k ~ 1 x Qk+l' k ~ 1 ' k ■■ V-« e Q^" 1 ^ 1 }. 

(2.11) 

We equip this space with surface constraints to obtain the full definition of the 
approximation space on the infinite pyramid: 

Definition 2.13. The kth order approximation space for H w (div, Qqo) is 



W (2) ' fe (f2oo) = {u E W( 2 )< fe | ri 2 ^ e P k -^[x,z] similarly on S 4 , 0oc ,z = 2,3,4}. (2.12) 

Again, we can identify the 2-forms in U^ 2 '' k (Cl 00 ') with vanishing normal traces on 
the faces of fioo. We denote this set by Uq (f2oo)- In the Appendix, we have written 
down a basis for U^ > ' k (p. 00 ). l and augmented it with shape functions for the faces. 

Since we want the divergence operator to be surjective as a map from W( 2 )' fc (fioo) 
to the associated approximation space of Lj^fioo), the approximation space for L 2 (i7oo) 
(considered as the space of 3-forms) consists of precisely divW^ 2 - ),,c (51 00 ). There is no 
longer any need to define an underlying space. 

Definition 2.14. We define the approximation space U^' k [iloo) f or -^^(^oo) 

as 

W (3) ' fe (ftoo) = Qi^'*" 1 '*" 1 . (2.13) 



3. The approximation spaces I4^' k (il) on the finite pyramid. We are now 

readily able to demonstrate the approximation spaces for the de Rham sequence on 
the finite pyramid. The construction is based on the approximation spaces on the 
infinite pyramid Cloo'. 



Definition 3.1. Let Q be the finite reference pyramid as defined in (1.1). Then, 
the kth order conforming subspaces on the finite pyramid Q are 

W (s) ' fc (!]) ^{(f^iueW 1 ^)}, s = 0,l,2,3. (3.1) 
Theorem 3.2. Let k be a positive integer. The finite dimensional spaces defined 



in (3.1 1 satisfy: 

W^ffljci? 1 ^), W (1) ' fc (f>) C if (curl, fi), (3.2) 

u^' k {n) c fr(div, n), w (3) - fc (f>) c L 2 (n). (3.3) 
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Proof. The proof follows from the definitions and properties of U^ s '' k (Qoo), the 
pull-back map <p, and Lemma |l.3| □ 

In the following subsections, we shall establish several useful properties of these 
spaces. The analysis will typically be performed for the approximation spaces on the 
infinite pyramid, where the basis functions are tensorial in nature, and hexahedral 
symmetries can be used, which allows for simple calculations in many cases. The 
properties of the pull-back operator will allow us to demonstrate the results on the 
finite pyramid. 



3.1. _ff 1 (il)-conforming approximation spaces. In this section, we demon- 
strate that the grad operator is injective on Uq° : (0) , the set of bubble functions on 
the pyramid. 

Lemma 3.3. Let U^' k {Vl) be the subset ofU^' k (Q,), consisting of functions 
whose trace onto the faces and edges of VL are zero. IfVv = for some v € U^ 0) - k (n), 
v = on ft. 

Proof. This follows from the divergence theorem. □ 

We can easily see that ^ 0),fc (O) = {((/T 1 )*?; : u G U^^oo)} ■ From the re- 
marks following p^] ), it follows that dim^ 0),fe (fi) = dim^ 0) ' fe (fi oo ) = (k - l) 3 . 
Note that from the definition o{l4Q ,k (Q 00 ) and the discussion in Section 1.3 the face 
traces of functions in ti^' k (Q) are compatible with those of neighbouring tetrahedral 
and hexahedral elements. Finally, the shape functions in the Appendix show that 
the edge traces are well-defined, and that edge traces can be specified in consistent 
manner. 



3.2. H (curl, O)- conforming approximation spaces. We shall establish that 
the grad operator maps U^' k (VL) into Z^( 1 )' fc (J7). This is an important step towards 
showing exactness of the diagram in |1.2| We then show that the curl operator is 
injective on a certain subspace of W^' fc (f2), which will be used in establishing uni- 
solvency of the edge elements on the pyramid. We will finally demonstrate a discrete 
Helmholtz decom posit ion. Note that from the definition of Z//g fe (^oo) and the dis- 
cussion in Section 1.3 the face traces of functions in U^ 1 '^^) are compatible with 



those of neighbouring tetrahedral and hexahedral elements. 



Lemma 3.4. The gradient operator is well defined as a map from U(°^ k (Q) into 



WW. fc (0). 

Proof. It is easier to work on the infinite pyramid. Recall that a basis for 
WW' fc (fi(x)) is given by functions of the form u a ^,c — fx+^n > where a, b and c are 
integers and a G [0, k], b G [0, k] and c G [0, k — 1] or u a ,b, c = yj^P • ^ e wm show 

that the gradients of each of these functions lie in W( 1 )' fc (ri). The result is trivial for 
c = 0. For c > 1, 



1 



l>a,b,c — 



(l + z) k 



fc+1 



ax a - 1 y b (z c+1 +z c ) 
bx a y b - 1 (z c+1 + z c ) 
K x a y b ((c~ k)z c + cz - 1 ) 
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If c < k — 2 then Vu Q ,6 jC G Qfc + i :fe:fc 1 x Qfc+i 1 x Qk'+i 2 - I R tne case c = k 
we can let r = a; a ?/ b in (2.7) and then the remainder 




VH 0lb)C - | r v * | = t^^t | teV" 1 **- 1 | , (3.4) 

which is in Q fe 7j' ' 1 x Q k \ j 1,-1 x <5^^ fc 2 . Finally, if c = fc then choosing r = 1 



in (2.7| suffices. Now use the definition of £/W> fe (Q) in terms of the inverse pull-back 
of functions in W( s )^(S7 oc ), and the commutativity of the grad with the pull-backs, to 
conclude the result. □ 

Note that the previous result also follows immediately from the (unproven) equiv- 



alent characterisations of the underlying spaces, ( |2.5| and (2.8). 

We need to identify the functions in U^' k {£l) with vanishing tangential traces. 

Definition 3.5. Define U^ ,k (Q) to be the subspace of functions in U (1 ^ k (£l) 
with zero tangential component on the boundary ofQ. Define ^ocuri(^) ^o 1 ^^) 
to be the set of gradient-free functions in IA^' ] ' k (fi) . 

From Lemma 1.4 we know that if u G U^' k {£l) has zero tangential traces on a 



particular face or edge of fl, then its pullback to f2oo will have zero tangential traces 
on the associated face or edge. This allows us to characterize U^' k (£l) and f ocuri(^)- 

LEMMA 3.6. Functions inU^' k {£l) can be represented as (0 -1 )* (u), where u G 
lA^' k ^(Oqo) have the form 

y(l-y)zqi \ z k-i /r x z\ 

x(l-x)zq 2 + j— k \r y z\, (3.5) 
K x(l-x)y(l-y)q 3 J V + z > \-r J 

i — s-\k—\,k—2,k — 2 s-\k — 2.k—l,k—2 s-\k — 2,k — 2.k—2 j /i \ /i \ 

where q G Q k+ { *Q k+1 *Q k+ { and r = x[l — x)y{l — y)p, p G 

Proof. It is easily verified that the functions u above have zero tangential traces 
on the edges and faces of £loo, and therefore their inverse pullbacks (</> -1 )*(u) belong 
to U^ ] ' k {VL). Note that it follows that dimU^ ,k {Q) = dimU^* (Q^) = k(k - l) 2 + 
k(k - l) 2 + (k— l) 3 + {k - If = 3fc(fc - l) 2 . □ 

The curl operator has a non-empty null space in IAq (fi), consisting of gradients. 
We can precisely characterize the complement of the gradients in Uq k (Q). 

Lemma 3.7. LetU^ := {v\v = (f'J'u^G w£url( fi oc)} where U^.^) C 
lA^' k (Slca) consists of functions u of the form 

y{l-y)zqi \ 

x(l - x)zq 2 , (3.6) 
K x{l - x)y(l - y)p) 

, _ ^,fe-l,fe-2,fe-2 _ ^,fe-2,fe-l.fe-2 ^,fe-2,fc-2r 1 rpi 7/( 1 ). fc /'r>\ r- 

where q x G Q k+ { ,q 2 G Q k+1 ,p G Q k+1 [x,y\. Then U^J. ml (il) C 

U^ ,k (Q), and the curl operator is infective onU^^ifl). In other words, gr&dUQ°^' k (Tl)P\ 
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Proof. The set inclusion j^ rl (fi) C (f2) follows by the definitions of 

^Oairi(^) an< ^ k {^)- To see that the curl operator is injective on Wg 1 ^^), we 

first show that the curl operator is injective on Wq 2^r(^oo)- The argument proceeds 
by contradiction. 

If k = 1 then WqJ^^Ooo) * s em Pty- Assume k > 2 and let u £ ^octjtiC^ 00 ) ^ e as m 



(3.6). Let either p or q 2 not equal to zero and write p — mffpri » r e 2 ' fc 2 ( a; : 2/)- 
Suppose that V x it = 0. From the cc-component, we obtain 

1 d d 

There is no z-dependence in r so we can factorise q 2 — f(z)g(x, y), where / e P k ~ 2 (z) 
satisfies 

d zf(z) _ 1 



dz{l + z) k + l (1 + 

This is impossible, and so p = q 2 = 0. A similar consideration of the y-component 
shows that q± = 0. We have just established that the curl operator is injective on 

^0 curi(^oo) - Since the pullback and curl commute, the curl is injective on curi(^)' 
□ ' 

We can now state a discrete Helmholtz decomposition for Uq (fi): 
Theorem 3.8. The discrete approximation space U£ ] ' k (ti) C ff(curl,n) of func- 
tions with vanishing tangential traces on dQ admits a Helmholtz decomposition. That 
is, if v € U<p ,k we can write v = + w with q G Wp ^' fc (f2) and w € ^o^uril^)- 

Proof. If q € U^' k (Q), it has zero trace on all the faces and edges of f2. Therefore, 
the tangential components of Vg are also zero on the fac es an d edges. We already 
know that grad maps U^> k (Q) into liW' k (Cl) from Lemma 



3.4 



and so it is clear that 



grad maps W (0) ' fc (fi) into U^ ] ' k (n). Injectivity of this map follows from Lemma 



Now we count dimensions. From Section 



and from Lemma 3.7 



3.1 



3.3 



we saw that dim^ 0) ' fe (O) = (k - l) 3 



dimW$£(fi) = dimi^ 1 ^^ = (2*c + l)(*c- l) 2 . 

From the same lemma, we know gradWg°^' fe (f2) H Wq^^O) = 0. Both of these are 
subspaces of ' fc (Q). So, 



dim {grad^ 0) ' fe (r!) U U^ lX (O) } = (2fc + 1)(* - l) 2 + (fc - l) 3 = 3fc(fc 



D 2 



which is the dimension of W^'^fi). Hence U^ ),k {n) = gradz4 0), *(fl) ©W^curl( Q )- D 
3.3. i?(div, fi)-conforming approximation spaces . In this section we shall 
establish that curlWW' fe (f2) C U^ 2 ^ k (fl). We then show that the div operator is injec- 
tive on a certain subspace of U^' k (fl), which will be used in establishing unisolvency 
of the divergence-conforming elements on the pyramid. We finally demonstrate a 
decomposition of this discrete space. 

Lemma 3.9. The curl operator maps elements of U^' k (fl) into I4( 2 > ,k ({l). The 
proof of this lemma is a calculation similar to the one in Lemma |3.4[ and is omitted 
here. 
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We now need to identify elements of U^' k (£l) which have vanishing norm al tra ces 
on the faces of the finite pyramid. Denote these by U^' k {Vt). From Lemma 



1.4 



we 

know that if T 2 n u = for some u g U^' k (Q) : then the pull-back L 2 n ^0*u = on 
the associated face of floo. This allows us to characterize U^ ,k (Vt) easily. 

Lemma 3.10. Functions inli^ ,k {VL) can be represented as where 
u G U^ ,k (yioo) have the form 

2t \ (x{l-x)xi 

2s I + I V(1-V)X2 | , (3.7) 

?X3 



(l + z) fe + 2 \ {1 + z){Sy+tx) 



where s = y(l - y)a, t = x(l - x)r, with \i G Qk+2 2 ' X2 G Qfe+2' fe 2 ' 2 > 

X3 G <5fc+2' fc ~ 1,fe ~ 2 ' °" G Q k ~ 1,k ~ 2 {x,y) and r € Q k ^ 2 ' k ~ 1 (x,y) 
Proof. 

It is easily verified that functions of the form (3.7| have vanishing normal com- 
ponents on the faces S^q^ of the infinite pyramid; their (inverse) pullbacks to the 
finite pyramid will thus have vanishing normal components on the faces Si si of fl. □ 
We note also that 

dimZ^'^fi) = dimWo^^ffioo) 

= k(k - l) 2 + k(k - l) 2 + k 2 {k - 1) + k(k - 1) + k(k - 1) 
= 3fc 3 -3fc 2 . 

We now present a subspace of U^ ,k {£l) on which the divergence operator will be 
injective. 

Definition 3.11. Define W (2 j-*(ft) := = (^- 1 )*(u),tt e W^v( n <»)} w/iere 



z fe-i / r y + 2t 



^o 2 div("oo) :=span{ - Z | r„ + 2a | } © span{ | )} (3.. 

^(l + z)(r xy + s v + t x )j 



and where r(x,y) = x(l — x)y(l — y)p(x,y),p <E Q k 2,fe 2 ,( = i(l-i)t,tGi 51 2 ( a; )> 
a = 2/(1 - S € P k - 2 {y), and X3 g Qj^a'* -1 '* -2 . 

LEMMA 3.12. T/ie divergence operator is injective on Uq 2 ^ ( fl) . 

Proof. We shall show that the divergence operator is injective on U^ k (fioo), and 
the desired result on f2 will follow by the properties of the pullback operator </> and 
the commutativity of </> with the divergence. Let u be as in (3.8). If V • u — 0, then 



(fc — l)z fc ~ 2 d 

= V ' u = (i + z )k+2 ( r *y +** + s ») + ai^ 3 )- 



We can factorize %3 = X^=o (i+fp+g gifo £/)■ We now compare coefficients of like 
powers of z. Since r, t and s have no dependence on z, we obtain 



= 



(k-l)z k 2 (r xy +t x + s y ) t d ( sr- z i+1 qi{x,y) 



(l + z) fc + 2 dz\^ (l + z) fe + 2 

Ji-1( • / k ~ 2 



(k-l)z k - 2 (r xy + t x + s y ) r^A z l+1 {i-k-l) + {l + i)z l 



(l + z) fe + 2 \^ (1 + 2!)* 
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This is impossible, unless 



(r xy + tx + s y ) = 0, qi(x,y) = 0. 



However, t only depends on x, and s only depends on y. From the form of r, it must 
be that r = = t = s. Therefore, V« 7^ for any non-zero u e Uq$^{£Loo)- Using 



the properties of the pullback operator, V ■ v = : 



for all v e U, 



(2),k 
O.div 



(Q). 



Finally, note that dimU^ 2 ^ (il) 
- 1. □ 



dimW 0div(^oo) - (k-l) 2 +2(k-l)+k 2 (k-l) 



'0,div 



fc 3 - 

Just as in the previous section, we can use Lemma (3.12) to exhibit a convenient 
decomposition of the discrete approximation space. 

Lemma 3.13. Any v € U^ 2),k (n) can be decomposed as v — V x w\ + W2 with 

Proof. Lemma 



3.9 



tells us that the curl operator maps W( 1 )^ fe (f2) into U^' k (Sl). 
Observe that if the tangential components of v are zero on some surface then the 
component of V x v that is normal to the surface will a lso be zero and so the curl 

we know that this mapping 



operator maps U Q j 
is injective. 



(fi) into Uq '' k (Q). By Lemma 



3.7 



By construction, 14$, j-^(fi) is a subset of U^'^in) and by Lemma 3.12 V • w 7^ 



(2),fc, 



for all w € ^Q 2 div(^)- Hence ^^y(fi)^^} 1 ( ]^ 1 (fi) is empty. We now count dimensions. 



We established in the proof of Lemma (3.12 ) that U^^Sl) has dimension A; 3 — 1 and 
from the previous section we know ^ocuri(^) nas dimension 2fc 3 — 3fc 2 + 1. Thus, 



dim(curlwW^(fi)UWg v fc (n)) 



which shows that U^' k {iX) = curing "^(fi) ®U^^l(VL). This establishes the desired 
decomposition. □ 

3.4. L 2 (£l)-conforming approximation spaces. We first note that the dimen- 
sion of U^> k (tt) = dimW^'^oo) = dim(Q k k ~ 1 3 ' k - 1 ' k - 1 ) = fc 3 . It is a straightfor- 
ward matter to determine that the divergence operator is well defined as a map from 
W( 2 ). fe (0) to W (3) - fe (ft). We record the result here in a lemma. 

Lemma 3.14. The divergence operator maps elements ofU( 2 ^ k (fl) intoU^' k (Q). 

Lemma 3.15. Any element u € U^' k (Q) can be written uniquely as 
u = V • w + A, we Ael. 



3fc 3 - 3k 2 



(2),fc, 



dimW (2) ' fc (fi), 



Proof. From Lemma 



3.14 



we know that divU^(Q) C U^' k (fl). We also know 



that the constants belong to U^ ,k (Jl). Now, dim(div Uq jj v (fl)) — k 3 — 1, which is one 
less than the dimension of U^ ,k (Q,). Now, suppose we could find w € div U^ k (fi) 

(2) k 

so that Vw = 1 on SI. By definition of ^- (fi), we know that w has zero normal 
components on the faces of SI. From the divergence theorem, this is impossible. Hence, 
we have shown that the constants are not contained in divW, 

divU. 



odiv(^)' an d therefore 



: £/( 3 )> fe (Q). This completes the proof. □ 
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4. First order elements on the pyramid. The approximation spaces l4^' k (Q) 
which we constructed include the elements presented by Gradinaru and Hiptmair [12"] 
as the special case k = 1. To demonstrate this, we will map the basis functions pre- 
sented in that paper onto the infinite pyramid, and demonstrate that these (pulled- 
back) elements belong to U^' k {p, 00 ). The properties of the pullback then allow us 
to conclude the set inclusions on the finite pyramid. The reason for this indirect 
approach is the tensorial nature of the approximation spaces on fi^ , which makes it 
easier to examine basis functions. 

• The lowest-order element: The basis functions for the 

element given in equation 3.2 of |12) are denoted 7Tj, i = 1..5. Set Tti — 



Tl = 
7T5 



(x-l)(y-l) 



1 + z 



7T2 = 



x(y- 1) 

1+Z ' 



7T3 = 



1 + z ' 



7T 4 



xy 

1 + z' 



z 

l + z 



It is clear that tt 1 E U^^iti^). 



• The lowest-order iJ(curl, fi) element: We proceed as in the case. 
Set 7i = 4>*"fi where the ji, i — 1..8 are the basis functions for the curl- 
conforming clement in [12 PI 



7i 



74 = 



77 = 



1 


(1 


+ z) 2 




1 


(1 


+ z) 2 




1 


(1 


+ z) 2 



1 



1 



(1 + z) 2 
1 

(1 + z) 2 

1 



73 = 



(1 + z) 2 



z(i-y) 

z(l — x) 
(l-y)(l- X ). 



76 = 




78 = 



(1 + z) 5 




(4.1) 



These are also the pullbacks of the basis functions for the first order curl 
conforming element given by Graglia et al. [T3]- Note that these are all edge 
shape functions. It is easy to see that ji are shape functions specified in the 
previous section for ff„,(curl, fioo) with k = 1. 

The lowest-order H(div) element: Set d = <fi*Ci, where Q, i = 1..5 are 

the divergence-conforming basis functions 



(1 + zf 



(1 + zf 




(2 = 



1 




c 4 - 



[1 + z? 



(1 + zf 




(4.2) 



3 There are minor typographical errors in 1 1 21 for two of the one-forms. Based on the preceding 
calculations in that paper, the correct expressions are 



76 



xz 

1-z 
xy , xyz 

1-z ~r (l-zp 



77 



\y 



1-z \ 

Z + 1-z 

xy , xyz 

1-z (1-z)'- 1 J 
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For completeness, we note that U^' 1 ^) consists of the constants, which map to 
multiples of nj^ya 011 ^ ne infinite pyramid. The above collections of functions are 



consistent with the definitions (2.4), (2.7), (2.10) and (2.131 



5. Conclusion. In this paper we have described high-order approximation spaces 
on a pyramidal domain. The construction is via an infinite reference pyramid which 
allows for simple, tensor-product descriptions of the approximation spaces. We en- 
sure compatibility across shared faces and edges with neighbouring tctrahedral and 
hexahedral elements by explicitly constraining the approximation spaces. We ended 
by demonstrating that the elements of Gradinaru and Hiptmair |12) can be recovered 
as the first order approximation spaces of our construction. 

In the companion article, we shall describe shape functions and degrees of freedom 
to complete the definition of the high-order elements. We shall also establish exactness 
of the discrete spaces. Since we are working with rational basis functions, it is not 
immediately apparent that high-degree polynomials are contained in these spaces; we 
shall demonstrate this polynomial approximation property as well in the companion 
article. 

Acknowledgement We gratefully acknowledge the contributions of Leszek Dem- 
kowicz, who suggested this problem. We would like to thank Leszek Demkowicz, Peter 
Monk and Paul Tupper for helpful discussions on the paper. The work of NN was 
supported by the Natural Sciences and Engineering Research Council of Canada, and 
the Canada Research Chairs program. JP was supported by a Natural Sciences and 
Engineering Research Council graduate fellowship. 

Appendix A. Shape functions on the infinite pyramid. 

We gather here, for reference, families of shape functions associated with the 
vertices, edges and faces for Z/A s '' fc (f2oo) for s = 0,1,2. We note that the shape 
functions presented below are not hierarchical, though a hierarchical construction is 
also possible. 

A.l. U^' k (n) shape functions. Since the approximation space W ( - ' ) ' fc (ri oo ) is 
invariant under the rotation, : ft^ — > O^, it is only necessary to demonstrate 
shape functions F a for a representative vertex, vertical edge, base edge and vertical 



face. Then, using (1.9) and the subsequent remarks, the inverse pullback of these to 
the finite pyramid will also be invariant under the rotation R. 
1. Vertex basis functions on Sloo 

( {l-x){l-y)} B - - - - ? k 

F Vl = < j- rv > , &iF V2 , F V3 , F Vi are similar and F V5 — 



(1 + z)fe j ,^ V2 ,, V3 ,, V4 _ —*„ B + 

There are a total of 5 vertex shape functions in U^' k (floo)- 

2. Basis functions for the representative vertical edge, e\\ 

f (1 - x)(\ - y)z a 1 

F e , = < ; r-= , 1 < a < k — 1>, and similar on the other vertical edges 

1 1 (Hz)' J 6 

There are a total of 4(fc — 1) shape functions associated with these edges. 

3. Basis functions for a representative base edge: 

h " =< (1 + z)* ' 



There are 4(fc — 1) such edge shape functions functions. 
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4. Basis functions for the representative vertical face S^q^: 

5. Basis functions for the base face, Bfi„: 

There are (fc — l) 2 such shape functions. 

A. 2. W^ 1 - ) ' fc (f2) shape functions. As in the previous case, we can perform this 
construction on the infinite pyramid and we will write G a — 4>*F a . 
1. Basis functions for the representative vertical edge, e\\ 



*(»-!) \ 
(a:-l)(y-l)(l + s)V 



(l + z) fc + 

2. Basis functions for the representative base edge: 



3. Basis functions for the representative vertical (triangular) face, Si.n^: 
We will use three types of basis function to construct F$ 1 : 



r , . /x tti (i + ^) 7i \ 

<W - S ,f , 1L > 7i,«i>0, 7 i+ai<fc-2 



(l + z)^ 1 









d a ,J (i - ri(i - ( rf; ) t + ' )> '°-'' fo i. o< M <*- 2 ; 



) 



G Sl =F Sl ,iUF Sl , 2 UF Sl 



It can be verified easily that G Sl C ZYW' fe . The size of the space spanned by 
Gsi is 

ifc(fc - 1) + \{k - l)(k - 2) + k - 1 = 2\k{k - 1) = dim^" 1 ) 2 - 

4. Basis functions for the base face, Bq^ 

( (V(l-V)x ai y*\ ) 

Gb = i jz , ,, , . , ai<fc-l,/3i<fc-2^ 



(l + 2) fc +! 









U | (1+ ^ )fc+1 \ x(l-x)x a *y^ J , a 2 <k-2,p 2 <k-l 

The size of Gb is equal to dim(2 fe_1 ' fe ~ 2 + dimQ k ~ 2 - k ~ 1 . 
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Lemma A.l. A basis for U^' k ((loo) * s given by the shape functions G a and 
functions of the form 

y{l-y)zqi \ (r x z^ 
x(l - x)zq 2 ) + {l + z)k+1 [r y z\, (A.l) 



,a;(l - x)y(l - y)q 3 



— r 



where q £ Q\. + \ k 2 ' k 2x Qk+i' k 2 x Qk+i' k 2 ' 2 and r = x(l — x)y(l — y)p, p E 
Q k - 2 > k - 2 [x,y]. 

A.3. Shape functions for U {2 ^ k (Vl). We will denote H a = 4>*F a . 

1. Basis functions for the representative vertical face, Si.n^: 

H Sl = { {l+ \ )k+2 (^(l-y)x a z b | , a,b>0, a + b<k-l 

The number of basis functions is k(k — l)/2. 

2. Basis functions for the base face, Bq^: 

"-^(j/ 0<a, b < k - iy 

Lemma A. 2. A basis for W^'^SIqo) is given by the shape functions H a and 
functions u G U^^^oo) which have the form 



It \ / x(l — x)xi\ 



{1 + Z) \{l + z){s y +t x )j \ z X3 J ( A ' 2 ) 
s = 2/(1 - y)o; t = x(l- x)t, 

withxi e Q fe+ 2 , X2 € Q fe+2 , X3 G Q fc+2 andcr e Q fc ^ 2 (a:,?/) 

and t e Q k ~ 2 ' k ~ 1 (x, y). 
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